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The scaling properties of a random walker subject to the global constraint that it needs to 
visit each site an even number of times are determined. Such walks are realized in the equilibrium 
state of one dimensional surfaces that are subject to dissociative dimer-type surface dynamics. 
Moreover, they can be mapped onto unconstrained random walks on a random surface, and the latter 
corresponds to a non-Hermitian random free fermion model which describes electron localization near 
a band edge. We show analytically that the dynamic exponent of this random walk is z = d + 2 
in spatial dimension d. This explains the anomalous roughness, with exponent a — 1/3, in one 
dimensional equilibrium surfaces with dissociative dimer-type dynamics. 
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I. INTRODUCTION 

Random walks provide the key to the scaling prop- 
erties of many physical phenomena. Some examples 
are: surface diffusion, wave packet spreading in quantum 
mechanics, equilibrium commensurate-incommensurate 
phase transitions in physisorbed monolayers on sub- 
strates, and one dimensional (ID) interfaces in equilib- 
rium and in stationary growing states In its ele- 
mentary formulation, a particle moves through d dimen- 
sional space by jumping during each discrete time step 
over a displacement vector x, according to a given (local- 
ized) probability distribution W(x). The fluctuations in 
the position, ri t — Ylt'=i &t'> a fter t time steps, obey the 
scaling form 



An t = (\n t - (ri t 
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and the spatial probability distribution P(ft, t) develops 
at large time scales into the Gaussian form 



P(n,t) = 
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■ exp 
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which is invariant under the scale transformation P — > 
b d P, t — > b z t and n — > bri. z is the so-called dynamic 
exponent, and has the simple value z = z rw = 2 in 
all dimensions, as is obvious from the linear character 
of the underlying diffusion equation. Scale invariance is 
generic to many other dynamic processes as well, but 
with typically non-trivial values for z. Kardar-Parisi- 
Zhang (KPZ) type surface growth [Q, and statistical pop- 
ulation dynamics like directed percolation and directed 
Ising type processes are examples of this ||. In such 
processes random walk (diffusion) arguments still play a 
central role, e.g., z rw typically appears within molecular 
field type approximations. 

The scaling properties of stochastic processes can be 
classified in so-called universality classes, according to 
the symmetries and conservation laws of the underlying 



dynamic rules. In analogy with equilibrium critical phe- 
nomena, dynamic critical exponents are universal quanti- 
ties, that are insensitive to varying details of the dynamic 
rule. For example, all random walks have z = z rw = 2, 
irrespective of the shape of the displacement distribution 
W(x). To break out of the z rw straight jacket, some- 
thing more drastic has to change. Examples of anoma- 
lous diffusion are: Levi flights, where the typical length of 
the displacement is divergent; correlated random walks, 
where the steps have long range temporal correlations || ; 
walks in a quenched randomness environment, like poly- 
mers in disordered media |^ ; and also diffusion on a one 
dimensional lattice, where the particles can move only in 
dimer form Q|. The latter leads to an infinite num- 
ber of conserved quantities, and density auto-correlation 
functions that decay as power-laws with anomalous ex- 
ponents. 

In this paper, we investigate the scaling properties of 
so-called even-visiting random walks (EVRW) on a d 
dimensional hypercubic lattice. During each time step 
(0 < t' < t) the random walker hops to one of its nearest 
neighbor sites with equal probability. However, unlike 
normal random walks, it is required to visit every site an 
even number of times before the walk terminates at time 
t. This gives rise to anomalous scaling. The even- visiting 
condition imposes a (mod 2 type) global constraint on the 
motion of the random walker, which it can satisfy only 
through correlated movements. The origin and nature 
of this type of anomalous scaling is completely different 
from those in the examples mentioned above. 

This study of EVRW's is complementary to our recent 
work on dissociative dimer-type surface growth in one di- 
mension (ID) ||. The surface grows and erodes by the 
deposition and evaporation of dimers only. Those dimers 
dissociate while on the surface (but do not diffuse) such 
that each monomer can arrive and depart with a differ- 
ent partner. This growth rule implies that the number of 
particles at each height level is globally (but not locally) 
conserved mod 2. Compared to conventional (monomer 



1 



deposition type) surface roughness, it imposes a global 
constraint on the fluctuations of the surface and leads to 
anomalous equilibrium roughness. The mod 2 conserva- 
tion of the particle number in dimer growth is equivalent 
to the even-visiting constraint in random walks, and the 
anomalous surface roughness is linked to the anomalous 
scaling of the EVRW. 

There exists a second completely different type of 
application of EVRW's. The time evolution operator 
of the random walk can be cast in the form of quan- 
tum mechanical non-interacting electrons moving in a 
random medium. The global EVRW constraint trans- 
lates into spatial and directional randomness of the hop- 
ping ampl itude s and a non-Hermitian random Hamilto- 
nian (Sec. IV E ). The spectral properties of such Hamil- 
tonians are a topic of growing research, in particular their 
localization-delocalization transition aspects [ ^o|Jlll . Our 
EVRW scaling involves only one electron and therefore 
relates to the scaling of the electron eigenstates near the 
bottom of the (almost empty) energy band. Those edge 
states have so-called Lifshitz tails, with essential singular- 
ities in the disorder-averaged density of states |l2[ . Our 
study of EVRWs elucidates the nature of the edge states. 

The EVRW problem naturally generalizes to Q-visiting 
random walks (QVRW) with a constraint that the num- 
ber of visits to each site must be multiple of Q. Diffu- 
sive motion of the QVRW describes the stationary state 
roughness of dissociative Q-mer growth, where a sur- 
face grows (erodes) by the deposition (evaporation) of 
a Q-mer. This conserves the number of particles at 
each height level modulo Q, which corresponds to the 
Q-visiting constraint. In our discussion we will focus 
mostly on EVRW and dissociative dimer type growth, 
but most results are easily extended to QVRW and the 
scaling properties turn out to be independent of Q. 

This paper is organized as follows. In Sec. ||, we 
review one dimensional dissociative dimer type dynam- 
ics Q , and present additional numerical results showing 
the anomalous roughness of the equilibrium surface. 

In Sec. Ill, we map this dimer growth model onto the 
ID EVRW and present numerical results for the anoma- 
lous diffusion in ID EVRWs. The finite size scaling type 
exact enumeration and Monte Carlo simulation results 
suggest that the dynamic exponent of the ID EVRW is 
equal to z ~ 3.0. We also devise an analytic scaling 
argument, a healing time argument, that predicts that 
the dynamic exponent is equal to z — d + 2 in d dimen- 
sions. The healing time argument presumes the existence 
of a crossover time scale Tf rcc ~ t 2 ^ d+2 \ within which 
the random walker is not influenced by the global even- 
visiting constraint imposed at time t. We numerically 
confirm that such a time scale is present in d = 1 and 
d = 2. 



In Sec. IV, we embed the EVRW into an Ising spin 
environment. A Sft = ±1 spin is assigned to every site. 
Each of them points up at time t = 0. Next, the particle 
moves as in a conventional unconstrained random walk, 
but the Ising spin at site ft flips with probability e (or 



not with probability / = 1 — e) each time the particle 
visits that site. The global EVRW constraint can now 
be represented by the requirement that all spins must 
be pointing up again at time t, i.e., by projecting out 
from all conventional RW configurations those that leave 
spins pointing down. The even-visiting walks condition 
obviously requires that e = 1, but it is advantageous to 
proceed with generalized values < e < 1. We call this 
generalized version the stochastic even-visiting random 
walk (SEVRW). Next, we construct the time evolution 
operator of the Ising spins and the random walker. The 
spin part is easily diagonalized and the result has the 
form of a non-Hermitian single particle quantum Hamil- 
tonian with quenched randomness. The latter can also 
be interpreted as the transfer matrix (thermal equilib- 
rium) of a polymer fluctuating in an environment with 
randomly placed defect-lines. 

In Sec. [y]we focus on one special point in the SEVRW 
model, e = f = 1/2. There the ID model is easily sol- 
uble. The dynamic exponent is exactly equal to z = 3. 
This point plays the role of stable fixed point in the sense 
of renormalization transformations in the EVRW phase 
diagram. The EVRW, although at the border, belongs 
to its basin of attraction. 

In Sec. VI we discuss the same issues as in Sec. 0, 



but now in the framework of the non-Hermitian quan- 
tum Hamiltonian. The anomalous dynamic exponent 
z = d + 2 relates to so-called Lifshitz tails in the density 
of states near the edge of the energy band. We co nclude 
with a brief summary and discussion, in Sec. VII. 



II. 



DISSOCIATIVE DIMER SURFACE 
DYNAMICS 



A. surface roughness 

Equilibrium and non-equilibrium properties of ID in- 
terfaces like crystal surfaces have attracted considerable 
interest jl3]| . Interfaces display intriguing scaling prop- 
erties and their growth dynamics is well understood in 
terms of a few universality classes. KPZ growth is one 
of the examples An interesting quantity associ- 

ated with interface roughness is the averaged width, W, 
i.e., the standard deviation of the interface height hi(t) 
(Z = !,...,£): 
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where L is the substrate size. The width satisfies the 
dynamic scaling relation 

W(L,t) = L a f(t/L^) , (4) 

where the scaling function f(x) behaves as 
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The stationary state roughness exponent a and dynamic 
exponent Zw — a/ (3 are universal quantities. 

In most growth models the structural properties 
of the depositing (evaporating) particles are explicitly 
or implicitly presumed to be geometrically featureless 
monomers. Nevertheless, the geometric features of the 
particle shapes can strongly affect the growth dynamics 
and the stationary state scaling properties || . 

B. dimer dynamics 

Consider a crystal built from atoms of type X. Assume 
that deposition and evaporation takes place in contact 
with a gas where the atoms only appear in molecular 
dimer form X2, and assume that such dimers can only 
land and depart from the interface when aligned with 
the surface. Upon deposition a dimer attaches to two 
nearest neighbor surface sites and loses its dimer char- 
acter. Upon evaporation, two nearest neighbor surface 
atoms form a dimer and depart from the crystal. This 
dissociative character of the dimers is the essential fea- 
ture leading to the anomalous surface roughness in the 
equilibrium state. 

We describe the ID surface configurations in terms of 
integer height variables subject to the so-called restricted 
solid-on-solid (RSOS) constraint, hi — hi +1 — 0, ±1. The 
dynamic rule is as follows. First, select at random a bond 
(1,1 + 1). If the two sites are not at the same height, no 
evaporation nor deposition takes place. If the two sites 
are at the same height, deposition of a dimer covering 
both sites is attempted with probability p, or evapora- 
tion of a dimer with probability q = 1 — p (see Fig. [l]) . 
Processes are rejected if they would result in a violation 
of the RSOS constraint. 

Surfaces growing according to such dissociative dimer 
dynamic rules behave fundamentally different from those 
following monomer-type growth rules. The latter, irre- 
spective of being in equilibrium or in a stationary growing 
state, display, with only a few very notable exceptions, 
the universal roughness exponent a = 1/2; as exempli- 
fied in the Edwards- Wilkinson (EW) @ and the KPZ |] 
universality classes. The universal value of a is under- 
stood from a random walk argument. To be precise, a ID 
surface can be mapped on the time trajectory of a parti- 
cle in ID by identifying the height hi at each site I with 
the particle position n t at time t = I. The steps in ID 
surfaces are uncorrelated beyond a definite correlation 
length. Therefore the particle performs a random walk 
with displacement fluctuations \n t — n#\ ~ (t — t') l f Zr ™ at 
large time scales. This yields the value of the stationary 
state roughness exponent a = l/z rw = 1/2. 

Dissociating dimer growth circumvents the random 
walk argument by means of a novel type of non-local 
topological constraint. The dimer aspect requires that 



the number of particles at every surface height level 
must be conserved modulo 2. The dissociative charac- 
ter of the dimers transforms this into a non-local global 
feature. This leads to various interesting phenomena. 
In equilibrium, the surface is rough but with anoma- 
lous scaling exponents ]9|Jl5[| . Out of equilibrium, while 
growing or evaporating, it always facets |9|. Moreover, 
when the model is extended by introducing a so-called 
reduced digging probability at flat segments, towards a 
directed Ising type roughening transition in the extreme 
no-digging limit, the roughness becomes even more com- 
plex ]9|,|l6|]. The non-equilibrium faceting aspects are al- 
ready well documented in Ref . || . Here we focus on the 
anomalous equilibrium roughness. 

C. anomalous equilibrium roughness 

At p = q the above dynamic rule satisfies the detailed 
balance condition and the stationary state distribution 
is a genuine Gibbs type equilibrium state. We study the 
dynamic scaling of the surface width via Monte Carlo 
(MC) simulations. The crystal size L is even, with pe- 
riodic boundary conditions, /i^+z = hi, and we use as 
initial condition a flat surface, hi — for all I. The 
surface width is measured and averaged over N s inde- 
pendent MC runs, ranging from N s = 5000 for L = 2 5 to 
N s = 500 for L = 2 10 . 

The results are shown in Figs. ^ (a) and (b). The sur- 
face width does not obey monomer growth type EW scal- 
ing with a — 1/2 and (3 — 1/4. The dimer surface width 
saturates slower ((3 < 1/4) and is definitely less rough 
in equilibrium (a < 1/2). Notice the large corrections 
to finite size scaling of the width in both the temporal 
and spatial domains. These prevent us from obtaining 
accurate values for the exponents a and f3 from simple 
log-log type plots of the width versus t and L. Instead, 
we define effective exponents 

a(L) = \n[W(mL, 00)/ W(L, 00)}/ \nm (6) 

and 

f3(t)=\n[W(L,mt)/W(L,t)]/\nm, (7) 

where m is arbitrary (we choose m = 2) and W(L, 00) 
denotes the saturated width. For a(L), we use data for 
L = 2 5 ,...,2 10 , and for (3(t), the data at L = 2 13 at 
times shorter than t < 10 5 where finite size effects are 
still invisible. The results are shown in Figs. || (c) and 
(d). We estimate 

a = 0.29(4), = 0.111(2). (8) 

and zw — 2.6(5), since zw = a/(3. The exponents are 
definitely different from those of ordinary equilibrium 
rough interfaces but the precise values remain uncertain. 

The mod 2 non-local conservation of particle number 
is clearly the most promising candidate for being the ori- 
gin of the anomalous scaling behavior; as confirmed in 
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the following sections. However there exist additional 
more local conserved quantities in the dimer dynamics. 
When a dimer desorbs or adsorbs, the surface heights 
at two nearest neighbor sites change by one unit simul- 
taneously. This implies conservation of the anti-Bragg, 
k = 7r, Fourier component of the surface height 

^=-7fE e_tt '^- ( 9 ) 
VL 1=1 

In other words, the dynamics is not ergodic; surface 
configurations with different values of hk— v are dynami- 
cally disconnected. Therefore the scaling properties may 
depend on the initial condition. Such types of effects 
are studied in Ref. jl5| in the context of dissociative k- 
mer growth in body-centered solid-on-solid type models, 
hi - hi+i = ±1. 



D. surface diffusion 

In our model the particles do not diffuse along the 
surface. In actual experimental settings, surface diffu- 
sion cannot be ignored. The k = it broken ergodicity 
is restored by diffusion, but the mod 2 conservation is 
preserved as long as diffusion across steps is forbidden. 
Such jumps to higher and lower levels are suppressed by 
so-called Schwoebel barriers j[7). This means that the 
anomalous surface roughness discussed here can be ob- 
served at time scales smaller than the characteristic time 
associated with jumps across steps, provided the other 
time scales are short (high surface deposition rates). 

To test the robustness of anomalous dimer roughness 
and to verify the essential role of the global mod 2 parti- 
cle conservation at each height level, we add to the dimer 
growth model diffusion of surface atoms within terraces. 
The surface is again described by integer height variables 
hi, subject to the RSOS constraint and periodic bound- 
ary conditions. The dynamic rule is as follows. Select at 
random a bond (I, and attempt with equal probabil- 
ity: a dimer deposition or evaporation just like above; or 
a monomer jump from site I to one of its nearest neigh- 
bor sites. The move is rejected if it would result in a 
violation of the RSOS constraint. Since the RSOS con- 
dition is imposed at every stage, jumps across steps are 
automatically forbidden. 

Starting from a flat surface at t = 0, the surface widths 
are measured for L = 2 5 , . . . , 2 9 . The results are shown 
in Fig. ||(a). They are qualitatively the same as in the 
absence of diffusion. The exponents a and f3 are deter- 
mined in the same way as in Eqs. iM) and (R), see Fig. || 
(b) : 

a = 0.31(3), = 0.115(5). (10) 

The finite size corrections to scaling are again very large. 
The exponents are slightly larger than in Eq. (||), but, 



within the current numerical accuracy we cannot distin- 
guish one from the other. 

We conclude that dissociative dimer equilibrium dy- 
namics represents a new universality class for interface 
roughness. Surface diffusion within terraces is irrelevant 
and this new universality class is characterized by the 
topological constraint caused by the mod 2 conservation 
of the number of particles at every height level. 



III. EVEN- VISITING RANDOM WALKS 
A. the model 

The above numerical study of dissociative dimer type 
dynamics clearly indicates that the equilibrium scaling 
properties of the interface belong to a different univer- 
sality class than conventional monomer type dynamics. 
We also identified the most likely origin of this: the con- 
straint that the number of particles at each height level 
must be preserved modulo 2 in a global non-local man- 
ner. The exact value of the exponent a is difficult to pin 
point from the MC results, due to strong corrections to 
scaling. To resolve this, we investigate in this section the 
properties of a random walk with the constraint that it 
needs to visit every site an even number of times before 
it terminates. This is the so-called even-visiting random 
walk (EVRW). 

Consider a random walker on a ID lattice, which is 
required to jump during each time step one site to the 
left or the right with equal probability, Uf+i = n t > ± 1. 
n t > denotes the position of the walker at time t' . The 
walker is demanded to visit every site n an even number 
of times after t time steps. 

We focus our presentation on the EVRW in one di- 
mension. The generalization to d > 1 is straightfor- 
ward and mentioned when appropriate. Moreover, it is 
natural to expand the EVRW into a Q-visiting random 
walk (QVRW) with the constraint that each site must 
be visited a multiple of Q-times. We obtained numerical 
results for Q > 3, but since we did not detect any differ- 
ences from the scaling behaviors at Q = 2 [Q, we limit 
this presentation to EVRW. 

The connection with dimer surface dynamics is self- 
evident. The probability distribution of EVRW repre- 
sents the equilibrium Gibbs distribution, i.e., the equi- 
librium state of a surface where all configurations that 
satisfy the mod 2 constraint are equally likely. There 
is one minor difference between our RSOS dimer model 
and the above EVRW. In the latter the particle is re- 
quired to make a hop during every time step, An = ±1, 
while in the RSOS dimer dynamics it is allowed to stay 
at the same site, An = 0, ±1. Figure [l| shows examples of 
both. This so-called body-centered solid-on-solid version 
of the EVRW is more compact and converges numerically 
faster. 
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B. exact enumerations 

The number of possible space-time configurations of a 
normal ID random walker is equal to Znwif) — 2*. The 
even- visiting constraint excludes most of those walks. 
It is of interest to know whether the total number of 
EVRW's still scales exponentially as Z(t) ~ /i 4 , and if 
so, whether [i remains equal to 2. For that purpose, we 
enumerate all EVRW's that start and return to the ori- 
gin (n = 0) after t time steps, using the exact (but not 
closed form) expressions, Eqs.(pd|) and ( |l2"| ) below, which 
were developed in Ref. Jl9[ in the following manner. 

Denote the number of steps to the right (left) from site 
n to n + 1 (n — 1) by r n (l n ). The number of visits of 
site n is equal to v n = r n + l n and the sum of all visits 
is equal to total number of time steps t = J2 n w « ■ The 
return-to-origin condition implies that r„ = l n +\, i.e., 
that v n = r n +r„_i = l n +i + l n , and that l n and r n must 
be even for all n, such that t = J>\ Vn ~ 2 J2 n Tn ^ s a 
multiple of 4 instead of 2. 

Every walk can be specified by the left boundary of the 
walk n m ; n , and m positive integer variables [s±, . . . , s m ]. 
The excursion m is defined by the distance between the 
right and left boundaries of the walk. The number of 
steps from site n to n + 1 is equal to r n = 2s„_„ min+ i 
with the understanding that s n < = for n' < and 
n' > m. The number of walks with the same set of posi- 
tive integers [s\, . . . , s m ] can be readily evaluated and is 
equal to B 
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t " 1 (2 Sj + 2 Sz+1 - 1)! 
2sx l\ (2s i+1 yp Si - 1)! [ ' 



for m > 2 and is equal to St si i — 2 for m = 1. The total 
number of the EVRW's is given by the sum 



t/i 



(12) 



=!{*i} 



where the prime in the second summation denotes the 
constraint that £ = 4^ . Sj, and the superscript in Z° 
represents the return-to-origin condition. 

Although analytically exact, this formula still involves 
infinite sums. Therefore we must resort to numerical enu- 
merations to determine the scaling properties. This has 
to be a finite size scaling type analysis because of the 
numerical upper limit for t. 

In Fig. [|(a), we plot Z°(t) as function of time for 
t < 140. The linear dependence in this semi-log plot in- 
dicates an exponential form Z°(t) ~ /i*. Next, we define 
an effective finite size exponent as 



Kt) = [Z°(t)/Z°(t-A)} 



(13) 



The corrections to scaling in Fig. |^(b) are strong, but a 
Neville type extrapolation analysis EQ] yields 



Despite the severe global constraint, the total number of 
EVRW's scales asymptotically in the same way as that 
of normal random walks with /i = 2. 

Figs. ||(a) and (b) indicate the presence of strong cor- 
rections to scaling. They are of an exponential form 

Z°(t) - 2 l e- atB . (15) 

as shown in Fig. ||(c). The slope yields 

9 = 0.34(2). (16) 

In Sec. 0, we will argue that the exponent 6 is a uni- 
versal quantity, and equal to the inverse of the dynamic 
exponent z of the EVRW; = l/z. 

We also performed an exact enumeration of the finite 
size scaling of the width of the EVRW. All configurations 
counted in Eq. (|ll]) have the same number of visits v n up 
to a constant shift in n. Hence, they all have the same 
width, W[{s t }} = (n 2 - n 2 ) 1/2 , with 



n = -J2(2i+l) Si 

i=l 
i=l 



The ensemble averaged surface width 

w2 = mhrr) E E w i{si}\ 2 s [su ... iSm] m 

K ' m=l { s ,} 

is evaluated numerically and plotted in Fig. ||(a). The 
roughness exponent a, W ~ L a , is estimated from an 
effective exponent 



Lf W{L) 
a{L) = 4 [w(L=T) - 1 



(18) 



see Fig. [3](b). Again, the convergence is slow, but the 
Neville type extrapolation yields 



a = 0.327(9) . 



(19) 



Within the numerical accuracy, this result is consistent 
with those in the two dimer growth models (with/without 
monomer diffusion) of the previous section, see Eqs. (|^) 
and @. 

The surface roughness exponent a is simply related to 
the dynamic exponent z of the EVRW as z = 1/a. So 
the above numerical result implies that 



z = 3.06(8) 



(20) 



fx = 2.000(2) 



(14) 



All the results of this section are checked numerically 
for Q = 3,4,5 in the QVRW model. We find no Q- 
dependence of the values of scaling exponents pa] . 
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C. Gaussian distributions 

We performed Monte Carlo simulations to determine 
the probability distribution P(n,t) for the EVRW, i.e., 
the probability to start at site n = and end after time t 
at site n. This was done by brute force. We simply gen- 
erated an ensemble of normal random walks and trashed 
the ones that did not satisfy the EVRW condition. The 
ratio decreases rapidly. For example, out of a total of 
2 x 10 9 normal random walks only about 600 walks sat- 
isfy the constraint at t = 500. 

The distribution function is shown in Fig. Q(a) and can 
be assumed to obey the scaling form 
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P(n,t) = —T{nlt^ z ) 



(21) 



with z the dynamic exponent. The best data collapse is 
obtained for 1/z = 0.32, as shown in Fig. ||. This value 
of z is consistent with the exact enumeration results of 
the previous subsection. It is also consistent with a di- 
rect evaluation of the second moment of the distribution 
function data, which yields that 



An 



5> 2 P(M)-(£ 

n \ n 



nP(n, t) 



1/2 



(22) 



scales as An ~ t 1 ^, with z ~ 3.3 as shown in Fig. 0(b). 

The functional form of the scaling function T is a sur- 
prise. It is of the form of T{u) — ^4e~ s l n l , as shown in 
Fig. |], with A ~ 1.98. This means that the probability 
distribution is Gaussian in nature, 



P(n,t) = Ar 1/z exp 



—B [n/t 



(23) 



This is surprising, because in other models with anoma- 
lous surface roughness, such as Levi flights, the probabil- 
ity distribution is certainly not Gaussian 

Gaussian distributions with z = 2 are characteristic 
for uncorrelated random processes. The appearance of a 
Gaussian shaped scaling function in the EVRW problem 
suggests us to search for an effective representation of the 
EVRW in which the correlation effects somehow trans- 
form away, with the possibility for an exact derivation of 
the EVRW dynamic exponent, possibly 2 = 3. This is 
the topic of the next section. 



IV. RANDOM WALKS COUPLED TO ISING 
SPINS 

A. defect spreading 

The even-visiting constraint is non-local in time. To 
keep track of this constraint in a local way, we can add 
an Ising field to a normal conventional random walk, i.e., 



a marker S n — ±1 to each site, that keeps track of the 
visits in the past. Initially at time t' = 0, all spins are 
prepared in the spin-up state. S n flips each time the 
random walker visits site n. The requirement that all 
spins are pointing up at t' = t, represents the EVRW 
constraint. The generalized distribution V({S}]n) t con- 
tains all the information we need, n is the location of the 
random walker at time t and {S} the spin configuration. 
P({+};ra)t is the EVRW distribution. 

Each down spin at intermediate time < t' < t repre- 
sents a defect, which needs to be healed at a later time. 
The defect area spreads in exactly the same way as the 
width of the conventional random walk, Wd ~ . We 
confirmed numerically that the defect distribution inside 
this cone is uniform in ID and 2D. This allows us to build 
the following healing time argument for the value of the 
EVRW dynamic exponent. 



B. defect healing time argument 

Divide the time interval t into two segments, Tf roe and 
Theal- For t' < Tf ree the random walker does not feel the 
constraint, diffuses freely, and leaves defects behind that 



1/2 
free " 



in 



are uniformly spread over a region of size £ ~ Tj 
order to satisfy the defect-free constraint at time t, the 
walker stops spreading and starts to heal defects during 
the second part of the walk, if > Tf roo . The typical dis- 
tance it needs to travel to heal a specific defect is of order 
£, and the time it takes the random walker to do that is 
of order t\ ~ £ 2 . The total number of defects Nd is of 
order £ d (d is the spatial dimension). Therefore, the heal- 
ing time scales as Th ca i ~ N d n ~ £ d+2 . Putting this all 
together yields a relation between the final time and the 
width £ of the EVRW. 



t = Tf r , 



T heal 



e 2 + e 



d+2 



(24) 



Thcai diverges faster than Tf ree , so we conclude that 

£^ f W+2) and z = d + 2. (25) 

The argument is more subtle in d > 2 due to the fact that 
the number of defects after time Tf rcc can not be larger 
than the total number of time steps, while the volume of 
the spreading cone, £ d ~ t^ 2 , diverges faster than that. 
This implies that in d > 2 the density of defects inside 
the spreading cone does not reach a constant. The num- 
ber of defects inside £ d is only proportional to Nd ~ £ 2 
instead of £ d . The time to heal one defect, n, however, 
changes as well, n is proportional to the time it takes to 
travel across the spreading cone £, times the probability 
to hit a defect while doing so, which is proportional to 
£ /Nd- The end result is that the healing time still scales 
the same as in d < 2, 



Theal ~ N d T\ 



Nd (£2jl)„£«H-2 

« Nd > 5 



(26) 
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We conclude that z = d + 2 in all dimensions. The value 
z = 3 in ID, is consistent with the numerical studies of 
the previous sections. This derivation is far from rigor- 
ous, but has the merit of being simpler than the ones in 
the following sections. 

The separation of t into two distinct time scales Tf rco 
and Thoai is artificial. Consider the average over all possi- 
ble starting positions of the random walker and all possi- 
ble spin configurations with periodic boundary conditions 
in the time direction (full trace). Then, the system be- 
comes translationally invariant in the time direction and 
the two distinct time domains should disappear. How- 
ever, Tf re e is still the natural crossover time scale in the 
problem. Consider the EVRW over time interval t. Mea- 
sure the width of the walk in a smaller time window t' 
inside t. For very small windows, t' <C t, the even-visiting 
constraint is invisible, and the width scales in the same 
manner as for a normal conventional random walk. This 
implies the following crossover scaling form for the width, 
£, of the EVRW 

£(t', t) = bt(b-*™t', b-'t) = t^G^'/Tfree) , (27) 

with b an arbitrary scale factor and z rw = 2. G is the 
crossover scaling function and Tf reo = t Zrw ^ z the crossover 
time scale. 

This crossover is important from a surface science per- 
spective. The time scale t corresponds to the charac- 
teristic length scale between impurities or other surface 
defects that act as effective lattice cutoffs. Depending 
on the experimental set up, like an X-ray beam width or 
STM scanning window that might be larger or smaller 
than this, one may measure the true asymptotic surface 
width scaling with a — l/(d + 2), or the unconstrained 
value a = 1/2. 

To illustrate the existence of this crossover time scale, 
we measure the spreading of the EVRW's 

A tnt , = [«T#» t -«r*»?] 1/2 • (28) 

((••■))t denotes the ensemble average over the walks that 
satisfy the even-visiting constraint at time t. Note that 
An in Eq. ( |22| ) is equal to A t n t . The spreading must obey 
the same type of crossover scaling form as in Eq. j27|), 

A t n t , = t l/(d+2)g( t , /t 2ttd+2)j ; (2Q) 

Monte Carlo simulations confirm this. We generate 
EVRW's over a given time interval t subject to the 
return-to-origin constraint, and record the time trajec- 
tories for < t' < t/2. Fig. ^| shows the spreading in (a) 
d = 1 and (b) d = 2. The crossover behavior is clearly 
visible in Figs. |^ (c) and (d). The data for different t col- 
lapse very well. As expected from Eq. (|2^), the scaling 
function Q(u) increases as u 1 ^ 2 in the short time region 
and saturates to a constant in the opposite limit. 



C. stochastic spin flip dynamics 

Consider a generalization of the EVRW in which the 
random walker flips the spin only probabilistically dur- 
ing each visit. The spin flips with probability e or is left 
unchanged with probability / = 1 — e. 

As in the deterministic EVRW problem, we require 
that at time t all spins return to the spin-up position. 
A more elegant and equivalent formulation of this is to 
require time-like periodic boundary conditions, because 
it suffices to demand that all spins at time t return to 
the same state as at time zero, irrespective of what that 
state might be, and the trace over all such initial condi- 
tions leads to periodic time-like boundary conditions. 

We call this model the stochastic even- visiting random 
walk (SEVRW). The deterministic EVRW corresponds to 
(e, /) = (1, 0) and the conventional RW to (e, /) = (0, 1). 

The purpose of this generalization is two-fold. On the 
one hand, it allow us to address the robustness of anoma- 
lous EVRW diffusion. On the other hand, and more im- 
portantly, there is an exactly solvable "decoupling point" , 
e = / = 1/2, where we can evaluate the anomalous dif- 
fusion scaling rigorously, see Sec. 0. 

D. non-Hermitian quenched randomness 

The master equation for the probability distribution 
reads 

V({S}; n) t+1 = | [V({S'}; n+l) t + V({S'}; n l) t ] 

+ t[p({sy,n +i) t +v({sy, n- i) t ] , (30) 

where configurations {S'} and {S} are related as S' n = 
—S n and S' m = S m for m =/= n. This can be cast in 
state vector notation, \V) t = J2{s} n n )t\{S}; n )> 

as \V)t+i = T \V)t with the time evolution operator 

f = i ^(e<7^ + /) [a) n a n+1 + a^-i] . (31) 

n 

The x-components of the Pauli spin operators, a x repre- 
sent the spin flips, and the fermion annihilation/creation 
operators a and a) represent the random walker. We 
have only one fermion in the energy band. 

The spin part of T is easily diagonalized since the a x 
do not couple to each other directly. Perform a rotation 
in spinor space to the eigenvectors, ^(|+) ± |— )), of &f l7 
and denote the eigenvalues as c n = ±1. In the rotated 
spinor basis, the operators a x become c-numbers, c n , and 
the time evolution operator reads 

T({c}) = ^^2{ec n + f)[a) n a n+1 + a],a„_i] . (32) 

n 

The initial all-spin up configuration becomes in the ro- 
tated spinor basis the linear superposition over all pos- 
sible {c n }. Each c n is either +1 or —1 at random and 
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does not evolve in time. The fermion (random walker) 
hops on a ID lattice with randomly placed defects, the 
c n = — 1 sites. The spin degrees of freedom transform 
into quenched random noise in the hopping probabilities. 
The time-periodic boundary conditions for the original 
spin variables translate into a quenched average over all 
defect configurations distributed uniformly. The wave 
function (probability distribution) is multiplied by a fac- 
tor (— e + f) each time the fermion visits a defect. Notice 
that the probability distribution can be negative when 
e > / for certain defect configurations. 

The generalization to Q-visiting random walks is 
straightforward. The eigenvalues become complex, c n = 
exp(i2wj/Q) with j — 1,...,Q, and (Q — 1) different 
kinds of defects appear with different random hopping 
probabilities. This type of generalization does not lead 
to any new scaling behavior of the probability distribu- 
tion of the random walker in the asymptotic limit |18|. 

The time evolution operator in Eq. ( p2|) resembles the 
Hamiltonian for an electron in a random medium. One 
fundamental difference is that T is non-Hermitian. The 
hopping probability from n to n + 1 is not Hermitian con- 
jugate to that from n + 1 to n. Non-Hermitian random 
Hamiltonians arise in various areas of physics. Stochastic 
processes, like random walks in disordered environments, 
have non-Hermitian time evolution operators. Equilib- 
rium systems with quenched disorder, like vortex line 
pinning in dirty superconductors fl(i| are described in 
the transfer matrix formulation by a non-Hermitian ran- 
dom Hamiltonian. Derealization transitions for such 
non-Hcrmitian types of disorder are different in nature 
from those in Hermitian systems, see e.g., Ref. 

This relation between non-Hermitian random Hamilto- 
nians and the EVRW is not new. It is presented typically 
starting from the non-Hermitian perspective. Our deriva- 
tion presented above using the reverse route has (in our 
opinion) the advanta ge o f being more transparent. To be 
precise, Cicuta et al recently considered a "roots of 
unity" model with Hamiltonian 



H = ^(& f „S„ + i + xilLi) 



(33) 



where b\ b* is a fermion operator and Xn = ±1 is the 
random variable with a uniform distribution. They re- 
late this non-Hermitian random Hamiltonian to the de- 
terministic EVRW. The disordered average of the trace 
of H l generates the EVRW configurations fll9|j . Alter- 
natively, the similarity transformation d n — C, n b n with 
|C„| 2 = 1 maps the time evolution operator in Eq. j3^ ) 
onto Eq. (d|) with c„ = Cn/Cn+i and Xn = c„c„_i. 



E. polymers in random media 



In the spin diagonalized form of Eq. (|32|), the single 
fermion is equivalent to a walker (fermion) in a quenched 
random environment. With probability 1/2 each site (n) 



is occupied by a defect, c n — —1, or not, c„ = +1. The 
probability distribution satisfies the recursion relation 

T({c},n) t+1 = \{ec n + f) 

[P({c},n+l) t +r({c},n-l) t ] (34) 

During each time step, V is multiplied with a factor ^ and 
with an additional factor f — e each time the walker lands 
on a defect site (recall that / + e = 1). This equation of 
motion does not preserve probability, and therefore we 
can not interpret it as a Master equation. The random 
walk nature of the problem is only restored after taking 
the quenched average over the c n randomness. 

Instead, we can interpret this equation of motion as 
the transfer matrix of a polymer wandering (but not back 
bending) on a 2D (n, i)-lattice with defect lines (at spe- 
cific n along the t direction). The partition function is 
equal to 



P(n,t) = 2 



-N, 



{c} walks 



2-V 



(35) 



with N s the number of lattice sites, v n the number of 
times the polymer visits site n in the specific walk under 
consideration, and ji = — log(/ — e) the energy associated 
with hitting a defect line. The prime in n' represents that 
we only sum inside the exponential over defect sites. 

The SEVRW interpolates between the normal ran- 
dom walk and the EVRW. At the random walk point, 
(e = 0, / = 1), the defects decouple from the polymer. 
At the EVRW point, (e = 1, / = 0), the summand in the 
partition function changes sign each time the polymer 
hits a defect line. 

Next, we can integrate out the defects altogether, be- 
cause the order of the two summations, the one over all 
polymer walks and the one over all possible defect line 
configurations, {c}, can be interchanged. (From the poly- 
mer perspective, the disorder is annealed, not quenched.) 
The trace over all defect configurations leaves us with 



p(m)= e 2-*n 



walks 



(36) 



with the product now running over all lattice sites m e 
N s . This leads us back into familiar territory. The 
SEVRW problem is now reformulated as a trace over nor- 
mal unconstrained RWs, but with Gibbs type weights 
giving each walk a different probability depending on the 
number of visits v n of every site. We could have started 
this way, because at = — 1 Eq. ( |3^ ) counts naturally 
only the EVRW, and at = 1 it counts all RW. For 
other values of \i the walks are weighted in a more com- 
plicated way, except at e = /, as we will discuss next. 
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V. THE EXACTLY SOLVABLE POINT 

A. reflective walls 

At point e = / the SEVRW is exactly solvable. Here 
the properties of the walk simplify in a different manner 
than at the EVRW point, e — 1, and the normal RW 
point, / = 1. The generating function representation of 
Eq.(|36|) reduces to 



P(n,t) 



walks 



-t-N v 



(37) 



with N v the number of distinct sites visited by that par- 
ticular random walk. The total number of walks is equal 
to 



(38) 



walks 



with h = log 2, and the summation running now over all 
walks irrespective of their end point. N v is also equal to 
the distance between the two extremal points reached by 
the RW. It is as if an energy is being assigned to each 
RW proportional to its space-time width. 

In the formulation of Eq. ( |35| ) the polymer is not al- 
lowed to cross defect lines (fi diverges), i.e., the problem 
factorizes in random sets of polymers on strips with finite 
widths. Similarly, the fermion time evolution operator 
reduces to 



T 



4 ^ 



+ 1) K^n+i + 



(39) 



boundary conditions. A heuristic evaluation of Z(£,t) 
runs as follows. 

For t < £ 2 , Z grows as Z ~ 2* just like normal random 
walks, but after this typical time scale the random walker 
begins to hit the boundary. It can only bounce back in- 
stead of having two possible futures (hopping directions). 
So compared to a walk in infinite space without reflective 
walls, the total number of walks is reduced by a definite 
factor each time the walker hits the wall. During time 
t, the random walker hits the boundary ~ i/£ 2 times on 
average. So one expects 



Z{^t)^2 t eM-at/e\ 



(41) 



with a is a constant of 0(1). The total number of con- 
figurations then scales as 



Z(t) 



J ^2 4 £e X p(-aV£ 2 -£log2) 



(42) 



The integral can be evaluated from the method of steep- 
est descent in the limit of large t : 



Z{t) ~ exp(-bt e ) 



(43) 



with 9 = 1/3 and b a constant. The maximum contribu- 
tion comes from £ m ~ i 1 / 3 and the power-law correction 
term follows in second order. 

The total number of walks returning to the origin, 
Z°(t), can be calculated in a similar way. The return- 
to-origin constraint reduces Z(£, t) by a factor of £. Wc 
obtain 



Z°(t) ~ 2H 1/6 exp(-fet 9 ) 



(44) 



The hopping probability to cross defect sites, c n = — 1, is 
zero. The defects act as hard core walls. The fermion is 
trapped and localized between two neighboring defects. 
These reflective walls are randomly distributed with a 
probability 1/2 to find one at every site without any spa- 
tial correlations. 

The probability to find in the quenched average the 
fermion within a blocked line segment of length £ is pro- 
portional to £2 - ^; because the probability to randomly 
place the fermion on a line segment of length £ is propor- 
tional to £, and the probability that such a line segment 
exists in the quenched average is proportional to 2~^. 
This allows us to calculate several quantities analytically 
in ID. 



B. total number of walks 

The total number of SEVRW walks can be reformu- 
lated as 



(40) 



where Z(£,t) is the number of possible normal random 
walks within a line segment of size £ with reflective 



C. spreading exponent 

The spreading, An(t), of the walker can be evaluated 
as well. First consider width u>(£, t) of a random walker 
trapped on a line segment of length £. Initially, for t < £ 2 , 
the random walker diffuses normally with w(£,t) ~ t 1 / 2 , 
until it realizes it is trapped. So w(£,t) saturates to 
£, and the spreading scales as w(£,t) = Cdi^ 1 /0 with 
g(x) ~ x for small x and g{x) constant for large x. The 
total spreading is the average of this : 



An(t) = (w(£,t)) 



f d£Z(£,t)£2-tw(£,t) 



/deZ(£,i)£2- 



(45) 



We use the method of steepest descent for large t, and 
again the maximum contribution comes from £ m ~ t 1 ^ 3 . 
This leads to 

An(t) - t 1 ' 3 , (46) 

i.e., z = 3 (since An(t) ~ t 1 ^), or after taking the 
crossover scaling into account, 

An{t) = i 1/3 <?(i 1/6 ) . (47) 

The crossover scaling dies out very slowly at large t, such 
that the corrections to scaling are large. 
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D. exponent identity 

We just established that the width scales as Art ~ t x l z 
with z = 3, see Eq. (E(j), and that the total number of 
walks has a correction factor exp [— bt e ] with 6 = 1/3, 
see Eqs. (|43| ) and (44)). We will demonstrate now that 
9 = l/z. 

The total number of constrained walks, Z(t), at the de- 
coupling point is given by Eq. (p8|). The average width 
of the random walk is equal to 



(N v [h}} 



d_ 

"dh 



\nZ[h] 



(48) 



Integrating this equation leads to the formal relation 



Z\h\ = Z[0] exp 



- / (N v [ti])dti 
Jo 



(49) 



with Z[Q] — 2*. It is reasonable to presume that (N v [h']) 
is continuous as function of h! . Then, according to the 
mean value theorem, the integral in the exponent is pro- 
portional to h{N v [h'o]) for < h'o < h. By setting 
h'o = log 2, we obtain 



Z(t) ~ 2 t exp[-a(V 1) )] 



(50) 



N v is equal to the excursion width of the walks, and pro- 
portional to An. Therefore 



Z(t) ~ 2* exp[-aAn(t)] 



(51) 



it remains surprising that the scaling properties of the de- 
terministic EVRW are the same as in the reflective wall 
model. 

The following intuitive derivation of Eq. (|l|) sheds 
some light on this. We expect that the total number of 
walks in every type of SEVRW is proportional to the total 
number of normal random walks 2* , times the probability 
that the Ising spin configuration satisfies the global con- 
straint. At the e — f point, the Ising spins flip randomly 
when their sites are visited. Therefore all spins inside 
the spreading cone are randomized completely and lack 
any spatial correlations. This means that the probabil- 
ity to find all Ising spins pointing up is proportional to 
exp[— aAn], which confirms Eq. (pTJ) . 

The extension of this argument to general SEVRW and 
the EVRW point in particular, requires that the distribu- 
tion of down spins is still uniform and that the spin-spin 
correlations are short-ranged in the large t limit. 

At the EVRW point, the random walker flips the spin 
at every visit. For large t, it is very likely that the num- 
ber of visits to every site inside the spreading region is 
even or odd with equal probability; we checked this nu- 
merically. Spin-spin correlations are the strongest at the 
EVRW point, but since this is a ID chain of Ising spins it 
is very unlikely that they can develop long range order of 
any type. We numerically measure the spin-spin corre- 
lation function, (S n S n + r ), and find exponential decay in 
the spatial direction; the correlation length saturates to 
a finite value for large t [Q . This explains why Eq. ( |5l| ) 
still holds at the EVRW point. 



Hence we conclude that the exponential factor in the par- 
tition function originates from the spreading of the walks 
and that 8 — l/z. 



E. universality 

Our numerical results for the EVRW model of the pre- 
vious sections agree with all the above exact results at 
the reflective wall point; see Eqs. (fl5|), (plf), (fl9|), and 
(pp|). This is actually somewhat surprising. 

It is relatively easy to argue that the scaling properties 
in the direct vicinity of the decoupling point e = / should 
be robust and universal, with the decoupling point act- 
ing as stable "fixed point" in the sense of renormalization 
transformations. At the decoupling point the fermion is 
deflected by the defects, while at e ^ f it can tunnel 
through them. This tunnelling is an exponentially small 
effect, see Eq. (|3(]). Passing through two defects is equiv- 
alent to passing through only one at a much smaller value 
of / — e, which means that under a rescaling of the spatial 
resolution the renormalized decreases towards zero. 

The normal random walk, at / = 1, and the deter- 
ministic EVRW, at e = 1 mark the natural horizons of 
the basin of attraction of this e = / fixed point. At 
these points, e _M becomes equal to ±1, respectively. So 



VI. LIFSHITZ TAILS IN RANDOM 
HAMILTONIANS 

A. density of states 

Let's return to the fermion time evolution operator 
Eq. (|32|), and examine the same scaling issues from that 
perspective. The number of walks Z (t) returning to the 
origin after t steps (n — 0) and satisfying the EVRW 
constraint can be written as 

Z o it) = 2 t 2 -NsJ2(0\(f({c})Y\0) 

{c} 

= 2* J dEp(E)E t , (52) 

where E is an eigenvalue of T and the disorder-averaged 
density of states is denoted by p{E). Since the operator is 
non-Hcrmitian, E is a complex number and the integra- 
tion runs over the complex E plane. Eigenstates near the 
band center are rather well documented for this type of 
non-Hermitian random Hamiltonians |jl^,^l|. However, 
we need to focus on the eigenstates near the band edge 
(at small wave numbers) since there is only one fermion 
in the system and our interests lie with its long time be- 
havior. 
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The nature of the eigenstates near the band edge, is 
rather well-known for Hermitian random systems. The 
density of states p(E) of these edge states exhibits an 
essential singularity, known as a Lifshitz tail |p^| . We 
review here an intuitive argument for the existence of 
Lifshitz tails and extend it to the non-Hermitian random 
SEVRW model. 



This exponential factor in the density of states near the 
band edge is known as a Lifshitz tail. Rigorous calcula- 
tions confirm its existence [Q. Moreover, the tails exist 
also in higher dimensions in the form of 

p(E) ~ exphalA^r^ 2 ] , (56) 
because, roughly speaking, pi then scales as . 



B. Lifshitz tails 

Consider a ID free fermion Hamiltonian with bond dis- 
order: 

H = — — ^ ^ r^m^n^m > (53) 
n.m 

where t n ^ m = t* m n are random hopping amplitudes be- 
tween sites n and m. This Hamiltonian is Hermitian. 
For simplicity, assume that the t n ^ m are nonzero only for 
pairs of nearest neighbor sites and take only the values 1 
and b (0 < b < 1) with equal probability. 

Without disorder, with all i„ jrra = 1, the energy 
band is trivial, E = — cos fc, with uniformly distributed 
wavenumbers, 5k = 2ir/L, in the range (— 7r < k < n). 
The states near the lower band edge, k ~ describe 
the large length scale behavior, and the density of states 
diverges as a power-law with the familiar van Hove sin- 
gularity 

p(E) ~ |AiT 1/2 , (54) 

in terms of AE = E — E e d ge ■ 

The eigenstates become localized in the presence of 
disorder. The probability to find a pure domain, i.e., a 
connected string of t n , m = +1, of size £ decreases ex- 
ponentially as 2~£. The crucial feature behind Lifshitz 
tails is that the states extending across the boundaries of 
pure domains do not contribute to the density of states 
near the edge, even in the presence of small tunneling 
probabilities (b > 0). In that case, the energy levels Ei 
in each segment are similar to those of a free particle in a 
box of size £, i.e., AE — E(k) — E e d ge ~ k 2 /2, with wave 
number spacing 5k — 27r/£; or, phrased in terms of the 
domain size £, |A£^| ~ (^/£) 2 for low-lying eigenstates 
with £=1,2,---. 

The distribution of first excited states p\ (E) between 
energy E and E + dE is proportional to the probabil- 
ity to find a pure domain segment with a size between 
£ and £ + d£, which is p 1 (E)dE ~ 2" s d£. Therefore, 
pi{E) ~ lA^l^expt-alAi;!- 1 / 2 ]. Similarly, for the 
l-th level, p e (E) ~ £\AE\- 3 / 2 exp[-a£\AE\-^ 2 }. The 
total density of states is the sum over all levels, but near 
the band edges, the contributions from higher levels yield 
only corrections to scaling. Hence the density of edge 
states is of the form 

p{E) ~ lA^r^cxphalA^r 1 / 2 ] . (55) 



C. Hermitian SEVRW model 

Let's now generalize this to negative hopping ampli- 
tudes. This may not be useful to real fermions in dis- 
ordered media, but is helpful to understand SEVRW's. 
Consider a Hermitian analogue of the SEVRW model 

T = - ^{ec n ^ n + f)aia m , (57) 

n,m 

where the sum is over nearest neighbor pairs and e = 1— / 
with < e, / < 1. The random variable c n ^ rn — c m , n can 
be either +1 or —1 with equal probability. So the hop- 
ping amplitude t n>m — ec„ !m + / is either +1 or (— e+f), 
and can be negative for e > /. The conventional Lifshitz 
tail argument applies to e < /. 

Similar to our earlier discussions, the c njm can be re- 
garded as eigenvalues of Ising-type spin flip operators 
<Jn,m- Unlike before, these Ising spins live on the bonds 
instead of the sites, e is the spin- flip probability when 
the walker (fermion) passes through the bond. Point 
(e, /) = (0, 1) corresponds to the normal RW model just 
like in the SEVRW model. However, there is an im- 
portant difference between the Hermitian and the non- 
Hermitian versions. The Hermitian formulation satisfies 
a self-duality relation between (e, /) and (/, e). The fol- 
lowing transformation on the creation/annihilation oper- 
ators 

n-l 

b n = a n Y[ c PiP+ i (58) 
P =i 

maps (e, /) onto (/, e). Therefore, the two limiting points 
(e, /) = (0, 1) and (1, 0) must correspond both to the nor- 
mal unconstrained RW. There is an even-visiting condi- 
tion at point (e, /) = (1,0), but it is imposed on the 
bonds. Unlike the site version, the bond constraint is 
automatically satisfied by all normal random walks re- 
turning to the origin. Consider a simple walk as exam- 
ple: walk ten steps to the left and then all the way back. 
When the RW turns around, it leaves a defect behind at 
the extremal point, in the site version but not in the bond 
version. On its way back it repairs all defects left behind 
during the first part of the journey, in both the site and 
bond versions. So in the bond version, all defects are au- 
tomatically repaired. In the non-Hermitian version (the 
original SEVRW model) the self-duality does not exist 
and (e, /) = (1, 0) is the anomalous EVRW problem. 
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At the decoupling point e = / = 1/2, the ID chain of 
Eq. (^) breaks up completely into randomly distributed 
finite segments, just like before in the non-Hermitian 
SEVRW. The hopping amplitudes t n ^ n are either +1 or 
0. These disconnected sections correspond to the pure 
domains in the Lifshitz argument and all states are com- 
pletely localized within those sections. The Lifshitz tail 
argument is exact at the decoupling point. 

The partition function, Eq. (|52|), is easily evaluated 
with the method of steepest descent as 

Z°(t) -2* J dE |A£|- 3 / 2 exp[~a\AE\-^ 2 - t\AE\] 

~ 2 t i 1 / 6 exp(-6t 1 / 3 ) . (59) 

As expected, we have exactly the same formula as in 
Eq. (Q) for the non-Hermitian decoupling point. Again, 
the dynamic exponent is equal to z = 3 in ID. 

In higher dimensions, the Lifshitz tails are of the form 

p{E) ~ exp[-a|A£T d/2 ] ■ (60) 

and therefore the partition function is proportional to 

Z°(t) ~ {2d) t eM-at d/(d+2) } . (61) 

Recall from Eq. (|oj) that 

Z(t) ~ Z°(t) ~ (2d)* exp[-a(AT„)] , (62) 

where (N v ) is the average number of distinct sites vis- 
ited by the constrained random walker after t time steps. 
Comparing these two equations yields 

(N v ) ~ t rf /( d + 2 ) , (63) 

and that every site is visited i 2 /( d + 2 ) times on average. 
This implies that (N v ) simply scales with the spreading 
volume (An) d ~ t d / z . Therefore, 

Z{t) ~ Z°(t) ~ (2d)* exp[-a< d / 2 ] , (64) 

with 

z=d+2. (65) 

This is the same result as obtained from the healing time 
argument for the EVRW, Eq. @. 

D. Lifshitz tails in the EVRW 

The Lifshitz tail argument also applies to the SEVRW 
time evolution operator, Eq. (|32|). Since T({c}) is non- 
Hermitian, the density of states is defined in the entire 
complex E plane. We focus here on the EVRW point 
(e, /) = (1, 0) where the distribution of states has a 
special symmetry property p9| . Apply the similarity 
transformation a n — e~™l 2 b n to the even sites and leave 
the odd sites invariant, a n = b n . T({c}) transforms to 



e -"/ 2 f ({/}) with c' n = c„ (-c„) for even (odd) n. Note 
that the disorder {c'} and {c} have the same distribution. 
Therefore, one obtains 

p(E) = p{e™l 2 E) . (66) 

This symmetry implies that there exist four Lifshitz tails, 
along the rays of arg(_E) = jir/2 with j — 0,1,2,3 at 
\E\ = 1. Each tail contributes equally to the parti- 
tion function Z°(t) apart from a phase factor exp(ij7ri/2) 
originating from the energy eigenvalue at each edge. So 
Z°(t) is equal to Eq. ( p9|) multiplied by the constant 
J3j=o exp(ij7rf/2). The latter is nonzero only when t is 
a multiple of 4, which is trivially true for EVRW's that 
return to the origin. We conclude that the dynamic ex- 
ponent for the non-Hermitian case is again z = 3 in ID 
and z = d + 2 in general dimensions. 

Finally, we can generalize to Q- visiting random walks. 
The analogue of Eq. ( |32"|) for d dimensional QVRWs is 
the time evolution operator 

1 ' 

T({c}) = — c n4«m (67) 

n m 

where n is a site of a d dimensional hypercubic lattice 
and the primed sum runs over nearest neighbor sites of 
given n. The random variable c„ takes equally likely the 
values exp(i2nj /Q) with j = 1, . . . , Q. 

The density of states has the symmetry property 
p{E) = p(e m /Q E), following the generalized similarity 
transformation d n = e~ m ^b n applied to one sublattice 
and leaving the others unchanged, a n = b n . Through this 
transformation, T({c}) picks up a phase factor e~ lrr ^ . 
There are 2Q Lifshitz tails, along the rays with aig(E) = 
jn/Q for j = 0, . . . , 2Q — 1 at \E\ = 1. Each tail con- 
tributes equally to the partition function Z° (t) except for 
the same type of trivial phase factors as in the EVRW 
(t is now a multiple of 2Q). The rest of the story is the 
same as for the EVRW, and the results are identical. 



VII. SUMMARY AND DISCUSSION 

In this paper we have investigated the scaling proper- 
ties of even-visiting random walks. The number of visits 
to each site by the random walker is required to be a mul- 
tiple of 2. This is a global constraint that leads to anoma- 
lous diffusive motion, of a novel type compared to more 
conventional ones such as Levi flights and correlated ran- 
dom walks. Using exact enumerations and Monte Carlo 
simulations, we find that the dynamic exponent is equal 
to z — 3 in ID. Surprisingly, the probability distribu- 
tion is not a stretched Gaussian (as for the other types 
of anomalous diffusion) but a simple Gaussian (with an 
anomalous value of z). We devise an healing time ar- 
gument which suggests that z = d + 2 in d dimensions. 
These results are verified numerically in ID and 2D. 
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We embed the even- visiting random walk into an Ising 
type environment, with an Ising spin at every site, where 
the random walker flips the spin at the site where it lands 
during each visit. Diagonalizing the spin sector of the 
master equation, translates the EVRW into a free fermion 
problem with quenched randomness. The time evolution 
operator takes the form of a non-Hermitian random-bond 
free fermion Hamiltonian. This leads naturally to the for- 
mulation of a generalization, stochastic even- visited ran- 
dom walks (SEVRW). The master equation for SEVRW 
can be reinterpreted as the partition function of a poly- 
mer fluctuating in an environment of randomly placed 
defect lines, and, after integrating out the randomness, 
as the equilibrium partition function of a polymer with 
an energy proportional to the width of the polymer con- 
figurations. 

The SEVRW model has a trivially exactly soluble 
point, the decoupling point where the polymer can not 
cross defect lines. At that point we can show rigorously 
that z — 3 in ID. Moreover, this point acts as a sta- 
ble fixed point in a renormalization transformation type 
sense in the SEVRW as a whole, such that z = 3 is 
valid in general. In the fermion interpretation, the same 
asymptotic anomalous diffusive properties of the EVRW 
determine the spectral properties of the non-Hermitian 
Hamiltonian near the band edge, in terms of so-called 
Lifshitz tails. This confirms that z — d + 2. 

The anomalous roughness we observed numerically in 
ID surfaces described by dissociative dimer-type dynam- 
ics was the starting point and motivation of this study. 
Such interfaces provide possible experimental realizations 
of EVRW's, like the roughness of steps on vicinal surfaces 
where the dynamics only allow attachment/detachment 
in the form of diatomic molecules. 

The scaling we found here is very robust. For exam- 
ple, they also apply to Q-mer type growth models. We 
established that the relation z = d + 2 remains valid for 
all values of Q in the Q-visiting random walk general- 
ization of EVRW's, and that the probability distribution 
still takes a Gaussian form. 

Random walks are a generic type of stochastic process. 
It will be very interesting to search for more novel types 
of scaling originating from RWs subject to global con- 
straints. 

During the final stages of preparing this manuscript, 
Bauer, Bernard, and Luck posted a preprint on the cond- 
mat archive pl| exploring the same type of connections 
between the EVRW and Lifshitz tails. Their results over- 
lap only partially with the research presented here. 
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FIG. 1. The upper panels show the deposition and evap- 
oration of a dimer. The lower panels show diffusion of a 
monomer. 
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FIG. 4. EVRW (solid line) and corresponding surface (dot- 
ted line) configurations. 




10 

10 1 \ 

10 

10" 



0.15 



(a). 



cn. 0.125- 




10 1 

10° 

0.5 
0.4 
0.3 
0.2 



(b) 



10' 



(d) 



10000 



0.035 
1 / L 



0.07 



100 



50 



[± - 



(a) 



2.0 



1.9 



(b) „ ( C ) 
= 10 - 



40 80 120 160 0.02 0.04 10 

/ 1// 



100 



FIG. 5. Exact enumeration of the total number of EVRW's. 
(a) lnZ° vs. t, (b) fi(t) vs. 1/t, and (c) log-log plot of 
-ln[2-*Z ] vs. t. 



FIG. 2. Monte Carlo results for the dimer model, (a) Time 
dependence of the surface width for L — 32, ... , 1024 from 
bottom to top. The straight line has slope 2/3ew = 1/2. 
(b) Saturated surface width. The straight line has slope 
2aEW = 1. (c) and (d) Effective values for /3 and a. 
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FIG. 6. Exact enumeration of the EVRW. (a) Surface 
width W. (b) Effective values of the roughness exponent a. 



10- 



0.035 
1 1 L 



0.07 



FIG. 3. (a) Time dependence of the surface width of the 
dimer model with diffusion, for L — 32, . . . , 512 from bottom 
to top. The straight line has slope 2(3ew = 1/2. (b) Effective 
values of a for with (□) and without (o) diffusion. 
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FIG. 7. Monte Carlo results for the ID EVRW. (a) Prob- 
ability distribution P(n,t) at t = 32(D), 64(o), 128(A), and 
256(v)- (b) Scaling of the spreading An. 
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FIG. 8. Scaling of P(n,t) for the ID EVRW at t = 32(D), 
64(o), 128(A), and 256(v) according to Eq. @. The best 
collapse is obtained with l/z = 0.32. The scaling function is 
assumed to be of the form J-(u) = Ae~ s '"' . A least square 
fitting yields A = 0.88, B = 0.62, and A = 1.98. 



FIG. 9. Crossover scaling in the EVRW. (a) and (b) 
Spreading of the EVRW in d = 1 and d = 2. (c) and (d) 
Scaling of the data according to Eq. (p9|). The broken lines 
have slope 1/2. 
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